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We study the polarization of the vacuum for a scalar field, {4> 2 ), on a asymptotically anti-de 
Sitter black hole geometry. The method we follow uses the WKB analytic expansion and point- 
splitting regularization, similarly to previous calculations in the asymptotically flat case. Following 
standard procedures, we write the Green function, regularize the initial divergent expression by 
point-splitting, renormalize it by subtracting geometrical counter-terms, and take the coincidence 
limit in the end. After explicitly demonstrating the cancellation of the divergences and the regularity 
of the Green function, we express the result as a sum of two parts. One is calculated analytically and 
the result expressed in terms of some generalized zeta-functions, which appear in the computation 
of functional determinants of Laplacians on Riemann spheres. We also describe some systematic 
methods to evaluate these functions numerically. Interestingly, the WKB approximation naturally 
organizes {4> 2 ) as a series in such zeta-functions. We demonstrate this explicitly up to next-to-leading 
order in the WKB expansion. The other term represents the 'remainder' of the WKB approximation 
and depends on the difference between an exact (numerical) expression and its WKB counterpart. 
This has to be dealt with by means of numerical approximation. The general results are specialized 
to the case of Schwarzschild-anti-de Sitter black hole geometries. The method is efficient enough 
to solve the semi- classical Einstein's equations taking into account the back-reaction from quantum 
fields on asymptotically anti-de Sitter black holes. 



I. INTRODUCTION 

Classically, matter influences gravity via its stress- 
energy tensor, T^ v , that appears as a source term in 
Einstein's equations. It seems natural to expect that, if 
quantum fluctuations of matter fields and the curvature 
are sufficiently small, this picture may remain valid also 
in the semi-classical theory, where quantum fields propa- 
gate on a given classical geometry. In this case, the quan- 
tum field would couple to gravity via its stress-energy 
tensor, and its 'back-reaction' effects on the background 
spacetime would be described by the semi-classical Ein- 
stein's equations: 

G^v — 87r(T M „) . 

Even though the precise range of validity of the above 
equations cannot be quantified without knowledge of the 
full quantum gravity, one may reasonably expect that, 
as long as back-reaction effects are locally small, semi- 
classical Einstein's equations may still provide an ade- 
quate description even if they lead to large cumulative 
effects, as in black hole evaporation. 

In general, calculating semi-classical back-reaction ef- 
fects is not an easy task, and it presents itself with a series 
of major difficulties. Ref. [l[ lists three of them. The first 
is related to the fact that in curved spacetime no renor- 
malization prescription can be uniquely given, and this 
generates an ambiguity in the expectation value of the 
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stress-energy tensor. Obviously, the semi-classical ap- 
proximation cannot resolve this issue. The second prob- 
lem is related to the fact that back-reaction effects intro- 
duce higher-derivative terms into the classical equations, 
possibly introducing spurious solutions. Therefore new 
criteria, beyond the semi-classical approximation, have 
to be invoked to select the physically relevant ones. A 
third problem stressed in [l[ is that: 'It is very difficult 
to compute (T^)'. Even for highly symmetric space- 
times, or for particularly simple vacuum states, it is a 
formidable task to obtain (T^) even by numerical ap- 
proximation. 

The problem of quantum back-reaction is particularly 
interesting in the case of black holes, and, more generi- 
cally, in spacetimes with horizons, because particle cre- 
ation may lead to significant effects that may cause the 
evaporation of black holes. For this reason, since the 
original discovery of Hawking 2] , the study of quantum 
effects from matter fields on black hole geometries has 
been the subject of much attention. 

What is really necessary for studying the effect of the 
back reaction is the full averaged energy-momentum ten- 
sor. However, (4> 2 ), which we will term for convenience 
vacuum polarization, already conveys much physical in- 
formation, with the bonus of a slightly less cumbersome 
computation. Candelas was the first to calculate the 
polarization of the vacuum for massless, minimally cou- 
pled scalar fields on a Schwarzschild geometry outside the 
horizon, and obtained a renormalized expression for ((j) 2 ) 
using the point-splitting method 3], The initial work of 
Ref. @ has been extended in various directions. Fur- 
ther analysis of the renormalized vacuum polarization 
has been done in [3| and the expectation value of the 



2 



energy-momentum tensor for the Hartle-Hawking state 
in the case of a conformally coupled scalar field on the 
same geometry has been obtained in 0, 0j- Other re- 
sults were reported in ,7j, where the expression for the 
vacuum polarization has been extended inside the black 
hole horizon. In Ref. Q Fawcett calculated numerically 
the energy-momentum tensor for conformal scalar fields 
and compared those results with Candelas' ones Q and 
with Page's analytical approximation 0, [l(|. Jensen 
and Ottewill have extended the above computation of 
the energy-momentum tensor to the electromagnetic field 

Since these initial works, the study of black hole radi- 
ance kept, and is still keeping, the interest in the prob- 
lem alive. Novel results have been obtained and extended 
the original ones to different background geometries and 
higher spin fields. Anderson considered the case of a 
massive field with general coupling and calculated the 
vacuum polarization on a Schwarzschild background fl2| 
and on an asymptotically-flat, static, spherically sym- 
metric spacetime [13[ , both in the case of Hartle-Hawking 
vacuum. These results have then been extended and 
the full energy-momentum tensor has been evaluated in 
Q3]. Based on the methods described in Ref. [3 111 El, 
Sushkov developed an analytical approximation and ob- 
tained (<fi 2 ) for Schwarzschild and wormholes spacctimes 
The analysis of quantum effects has also been ex- 
tended to consider various aspects of quantum effects in 
rotating black hole geometries in [I?], [H, Ell and re- 
cently further work in the case of Reissner-Nordstrom 
[20I ] and lukewarm black holes has also been studied [H[ • 

More recently, renewed interest in the quantum back 
reaction on black hole geometries has arisen in the 
context of the AdS/CFT correspondence. The five- 
dimensional brane world model proposed by Randall and 
Sundrum with an infinite warped extra-dimension [22j 
has been noticed to be possibly equivalent to the four- 
dimensional Einstein gravity with CFT corrections [23| . 
If we believe this correspondence, gravity in the brane 
world can be understood by inspecting the quantum 
back-reaction due to CFT. Applying the correspondence 
to black holes, a conjecture has been proposed: a 
brane localized black hole in Randall-Sundrum brane 
world model cannot stay static 0, l25| . In this con- 
text, the energy- momentum tensor due to a free scalar 
field on a four-dimensional Schwarzschild background 
in the Boulware state, which is compatible with the 
asymptotic flat condition, was computed, and it was 
confirmed that the quantum back-reaction diverges on 
the horizon (26|. Although there have been many at- 
tempts to examine whether this conjecture is correct or 
not [13, [H [M [HI, [H, [H [H, we have not reached 
any definite conclusion yet. The original conjecture is 
about asymptotically flat brane geometries. One of the 
present authors suggested that one may be able to shed 
new light on this issue by exten ding the setup to asymp- 
totically AdS brane geometries |35l |36|. What will hap- 
pen on the CFT side has already been discussed by Hawk- 



ing and Page many years ago [371 ]. In this case it is ex- 
pected that a static black hole solution exists even after 
taking into account the quantum back-reaction if the size 
of the black hole is sufficiently large. However, such a 
self-consistent quantum black hole solution has not been 
realized numerically so far. 

The above reasons provided enough motivation for us 
to consider the effects of quantum fields on asymptoti- 
cally anti-de Sitter (AdS) black hole geometries. We will 
not solve the full back-reaction problem in this notes, 
rather we will focus on the computation of the vacuum 
polarization for massive, non-minimally coupled scalar 
fields, leaving the details of the full energy-momentum 
tensor calculation and of the effects of the quantum back- 
reaction to a forthcoming publication [38| . 

The strategy presented in this paper is valid for a gen- 
eral spherically symmetric, asymptotically AdS space- 
time, but, for brevity, we focus our attention to the case 
of Schwarzschild-anti-de Sitter (Sch-AdS) black holes, de- 
scribed by the following form line element: 



ds 2 = -f(r)dt 2 + f-\r)dr 2 + 



r 2 dn 2 



where 



1 2M 

f(r) = l-— + r 2 k 2 , 



(1) 



(2) 



= G~ 1/2 is the Planck mass in natural units and k 
is the AdS curvature length. The horizon is located at 
r = ?7j, where f(rh) = 0. Customarily, the time coordi- 
nate is complexified, r = it, rendering the metric positive 
definite outside of the horizon. The apparent singularity 
at the horizon is removed by regarding the r coordinate 
as periodic with period (3: 



<][_ 

dr 



that, for Sch-AdS black holes, becomes 



= 4tt 



1 



Analogously to the asymptotically flat cases, asymptoti- 
cally AdS black holes have a characteristic temperature, 
T = /3 _1 , and an intrinsic entropy equal to a quarter 
of the area of the event horizon. However, unlike the 
asymptotically flat black holes, there exists a minimum 
temperature which occurs when its size is of the order of 
the characteristic radius of the AdS space. As a working 
example, we will consider the above AdS black holes and 
the numerical analysis will be specialized to this case. 

The results we will present are technical and mainly 
follow the direct procedure outlined by Candelas in the 
case of Schwarzschild black holes. The Green function is 
expressed as a sum of the products of the normal modes, 
therefore the first step in the computation requires the 
homogeneous solutions of the matter field equations to 
obtain a complete set of solutions. The choice of the 
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normal modes selects a specific vacuum state. Although 
the wave equation is separable in the spherically sym- 
metric case, it is not in general analytically soluble, and 
some numerical analysis is necessary. However, to prac- 
tically perform renormalization of the Green function, 
which originally diverges quartically, we need some ana- 
lytic approximation to the normal modes. For this pur- 
pose, the WKB approximation is useful. In the following 
we illustrate how to obtain the WKB approximants for 
the normal modes, highlighting the differences with the 
asymptotically flat case. 

After this initial part, we discuss the renormalization 
of the Green function. To subtract the divergent geomet- 
rical counter-terms, we first need to regularize the diver- 
gent sum. For this purpose, an appropriate procedure is 
the point splitting regularization, developed in general by 
Christensen [3^|. Relaxing the regularization after sub- 
tracting the divergent pieces, we obtain a renormalized 
expression for the coincidence limit of the Green function, 
which is formally divided into two parts: one amenable 
to analytical computation and one that requires numeri- 
cal evaluation. The latter part can be arranged so as to 
get fast convergence in the mode sum. 

The point splitting method guarantees the regularity 
of the Green function. Customarily the explicit proof 
of the regularity is not illustrated. Despite of the ap- 
parent superfluousness, we explicitly perform this com- 
putation, which, aside of being a non-trivial check of a 
complex algebraic calculation, will turn out to be useful 
to convince the reader of our systematic way of eval- 
uating the renormalized vacuum polarization. In fact, 
the renormalized propagator is made of various diver- 
gent parts, which, combined with each other, lead to a 
regular expression. Recognizing the structure of the di- 
vergences for each of these pieces allows us to perform 
a piece-by-piece renormalization. This will allow us to 
handle the numerical computation in a more convenient 
way. In fact, as we will see, the analytical part of the 
propagator can be expressed in terms of certain regu- 
larized (—functions, which appear in the computation 
of determinants of Laplacians on Riemann spheres [40| . 
The conclusive computation of the vacuum polarization 
is illustrated in Sec. fVl where the numerical and analyti- 
cal calculations are presented in detail. Our final remarks 
close the paper. 



II. GREEN FUNCTION 

The Euclidean Green function for a scalar field satisfies 
the following equation: 

(□ - m 2 - £R) G E (X, X') = -y/g8{X, X 1 ) . 

Choosing as vacuum state the Hartle-Hawking one, it can 
be expressed as 



G E (X,X') = 1 ^y^n(r-T')//3 



4tt 



which is characterized by the regularity and the pe- 
riodicity along the t direction with period (3. Here, 
cos 7 = cos 9 cos 9' + smd sm9' cos(0 — </>'). Writing the 
delta function as 



one can find that the function G n \ satisfies 

2 2 

1(1 + 1)- 



— U 2 f— 

dx \ dx 



f 



(m 2 + £R)x 2 }G n i{x, x') = -S(x - x 1 ) 



(4) 



In the above equation we have rescaled the radial coor- 
dinate, the mass and the curvature as: 



x = kr 
and defined 

with 



rh — m/k , 



R = k~ 2 R 



an 



(5) 
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which for Sch-AdS black holes takes the following form: 

+ 3x h 



2 \Xh 



where Xh = At/j. The function G n i(x,x') can be written 
in the usual way as the product of the two independent 
solutions of the homogeneous equation associated with 
Eq. Q: 



-{rh z + ^R)x z jtp(x) = 

Using tortoise coordinates, 

dx* = dx/ f , 
the homogeneous equation can be written as 
d 2 , {1(1 + 1) , /' 



(6) 



nr -\-'Jl j ./• (x<p)=0, 
(7) 

where ' represents derivative with respect to x. It is 
instructive to examine the behavior of the solutions near 
the horizon and at infinity. In the near horizon region 
the solutions are determined by 



dxi 



(xip) - 



(8) 
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This leads, as in the asymptotically flat case, to expo- 
nential solutions, 



<p ~ e 



±UJ n X* 



fx , 



(9) 



one of them being regular (at the horizon). At infinity 
(large- a;) one finds 



d 2 i d ! .., 

dx 2 x dx x 2 



^~0, (10) 



which, differently from the asymptotically flat case, ad- 
mits solutions of the form 



K-R) 



(11) 



Going back to the original equation Q, its solutions 
should be chosen by specifying their asymptotic behavior 
so that they are regular at infinity and on the horizon. 
We indicate as pf(r) the solution regular on the horizon 
and as qf(r) the one regular at infinity. The WKB ap- 
proximation of these solutions will be discussed in detail 
in the next section. Notice that, in the Schwarzschild 
case, Eq. (0} takes the form of Heun equation, with two 
regular singular points and one irregular singularity at 
infinity. The 4D Sch-AdS case is different. Eq. (g]) has 
five regular singular points, with infinity being a regular 
singularity. 

The n = 0, conformally coupled case is somewhat spe- 
cial, because at the horizon the zero order WKB approx- 
imant vanishes in the limit £ = 1/6 and m = 0. In the 
n = case, the behavior at the horizon can be easily 
understood by re-expressing the mode equation in terms 
of the the logarithmic derivative of ip near the horizon. 
It takes the form 



x 2 / 



, dhup 
dx 



1(1 + 1) - 2x 2 = 



from which it is easy to read off the behavior at x = x^: 
_ 1(1 + 1) - 24 



x h + Zx\ 



Standard procedure leads to the following expression 
for the Green function: 



G E (X,X') 



1 
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pY(x < )qf(x > ) 



Pi(cosj) 



x 2 f (q"dpf/dx — p"(x)dqj l /dx) 



where e = k(r — t'). The above expression is yet formal, 
because we have not found any explicit solution for the 
mode functions. In fact, as one can easily guess, such a 
solution can only be found numerically or by using some 
approximation method. 



III. WKB SOLUTIONS 

A convenient way to obtain an analytic expression for 
the solution of Eq. (J6j> is to use the WKB approximation. 
As we will see, this approximation is suitable and suffi- 
cient to perform all the renormalization procedure that 
we will need later. 

The WKB method can be implemented by writing the 
solution as 



/ \ — 3/2tia— n ± f x W(x')h(x')dx' 

ip(x) = x ' W v e Jx h v ' v ' 



(12) 



with 1] > 0. The + sign refers to the solution regular at 
the horizon, which we have called pf(x) and the — sign 
to qf(x). The overall factor x~ 3 / 2 is multiplied so that 
the asymptotic form at x — > 00 becomes compatible with 
the WKB ansatz: tp - x -| e ±lV9+4(m 2 +««) i°g^ Substi- 
tuting the previous ansatz in the homogeneous equation 
associated with (|4]), one easily finds 







w ± 



3 

xli 



h 2 



(x 2 f)' 

x 2 fh 



3 

xh 2 



(x 2 f)'\ W 



x 2 fh 2 J W 



W± 
r\ W" 



W 



?7(1 + 


v) w' 2 
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3(x 2 f)' 1(1 + 1) 


h 2 

f 2 h 2 


w 2 

, .2/ 

+ m —z - 

x z 


4a; 2 h 2 
- f 

X 


2x 3 fh 2 x 2 fh 2 



We choose the coefficient 77 and the function h so as to 
cancel the term with (zb)-signature. Hence, we have rj = 
1/2 and h — x/f, and the WKB equation simplifies to 



W 2 



ai 



w 
W 



a 2 



W 2 
W 2 



W" 

as W , (13) 



where 















3//' 


3f 2 


2 X 3 


4 a; 4 



^7, 



rh 2 - 



and wc have defined 



a 1 



J_( f _L 

2x 2 V x 



fl 2 = -T^J 



It 

4 x 2 



£i?4 

ar 



a-3 



t_ 

2x 2 



As it is well-known, the WKB method expresses the so- 
lution iteratively as 

W = W {0) + W {1) + W {2) + ■■■ . 

The leading order term is 

the next-to-leading order correction is computed by 
adding the derivative terms in (fT3^) evaluated for W = 



W<®, and so on, iteratively, to the desired order. In the 
following, everything is calculated up to next-to-leading 
order, i.e. W^K 

To be more precise, we introduce the following trun- 
cated WKB approximation of 1/W: 



1 1 

W n (z) = $(z)i/2 ~ 4$(z) 3 / 2 



(14) 



where the function ^(z) conveys the next-to- leading or- 
der WKB corrections and is defined by 



where 



a 2 



Q2 



$2 

$2" 



a 3 - 



a 2 - a 3 



5/ 2 



(15) 



(16) 



Going to higher order in the WKB expansion means 
adding higher order terms obtained by reiterating the 
above procedure. The explicit WKB expansions are 
lengthy and we will not report them here, but they can 
easily be handled by any symbolic algebra manipulation 
program. 

One can easily check the quality of the WKB approx- 
imation by comparing it with the numerical solutions. 
It turns out that the numerical solutions are well repro- 
duced by the WKB approximation already for small val- 
ues of n and I, as it is easy to argue. In fact, expanding 
W to next-to-leading order, one can easily see that the 
remainder in 1/W is of order 0(l~ 5 ,n~ 5 ). 



where 



W n {l) 



2x 

7 



G?ln<7™(x) dlnpV'(x) 



dx 



dx 



In the above expression the coincidence limit cannot be 
taken in a straightforward way because the sums over 
I and n are divergent. The divergence related to the I 
summation is not a serious one and can be by-passed 
easily by noticing that the Green function is a finite and 
well-defined object as long asl/ X 1 , i.e. e ^ 0. The 
trick commonly used is to add a multiple of 5(t — t') (and 
its derivatives if necessary) which, as long asr^r', does 
not alter the result. Basically, we can add terms of the 
form: 



k a \ - 



where a function Ri (x) independent of n can be chosen 
freely. By looking at the asymptotic behavior of the so- 
lution for large I, which can be obtained from the WKB 
approximation, 



Wn{l)~W n (l)~(l + l/2)f 



0(0 + 1/2)- 



we can easily find that choosing Ri(x) — \/{x\ff) is suf- 
ficient to remove the divergence in the summation over 
i, leading to 



IV. COINCIDENCE LIMIT 

Our main goal is to compute the coincidence limit of 
the Green function which will provide us the expression 
for (4> 2 ). Here we follow the same procedure that has 
been used in the Schwarzschild case. First take the par- 
tial coincidence limit, by setting r = r' and f2 = fi'. In 
this case the expression for the Green function simplifies 
to 



k 2 a °° 
G B (X,X') = -^Y. 



Anae 



-1/2) 



1=0 



x 



■ 3 W n (l) 



(17) 



I 2 00 

Oe(x,x') = ^J2^ 

n=0 



1=0 



1 + 1/2 



X 3 W n (l) Xy/f, 



(18) 



The above expression is still (ultraviolet) divergent in 
the coincidence limit due to the summation over n. In a 
general spherically symmetric space-time, the expression 
that needs to be subtracted from the Green function, 
before taking the e — * limit, is known to be given by [H, 

in 



Gdi 



div 



16ir< 



4 
^7 







~m 2 fe 2 ~ 




K 


R}j (in 




4 



1 

12/ 



lfl 
6 dr 2 



1_# 

3r dr 



(19) 



where is the Euler's constant. As shown in Ref. [13| . a more suitable way to combine the above counter-term with 
the original divergent expression for the Green function is to re-express (|19p . by using the Abel-Plana summation 
formula, as 



6 



where 



Ai 



(21) 



and A2 = A2,i + A2,2 + A2,3 with 

.,2 



TO 

2o7 



In (to 2 /) — In (a + •v/a 2 + TO 2 / ) H ~ 

v ' a V / I671" 2 

+7/ ^nry (V« 2 (i + it) 2 + ™ 2 f - \/« 2 (i - ^) 2 + ™ 2 /) , 



(g-l/6)i? 
2 

+ 2z 



1 



In 



^/a 2 + m 2 / " 
eft / 1 



(a + \/a 2 + to 2 /) + — In (to 2 /) 



1 \ yja 2 (l + it) 2 + rh 2 f ^ol 2 (\ ~ it) 2 +m 2 f 



A 23 = —(—j' 2 --f"--f'-m 2 

After subtracting (f2T)|) from HHJ), we obtain an expression for the renormalized Green function, and its full coincidence 
limit e — ► 0, which provides the renormalized vacuum polarization, can be taken: 



HX))=G% en) 



7 2 ( 00 / 1 

( X„Y,^|-A 1 -A J + g(( i + 1/2 ,- 5 ^ 



1 



(/) xy/J 



00 , 



w-„(0 aV7/ / 



' + + f m 2 + (f - l/6)i?) — 



.(22) 



In the above relation we have separated, for convenience, the contribution from n = 0. As shown in |3[, this n = 
contribution vanishes in the Schwarzschild case, but it does not in the asymptotically AdS case. 



VACUUM POLARIZATION 



A. Summation over I and Regularity of the Green function 



We explain how we perform the summation over I in (|22|) in this subsection. At the same time, we explicitly 
demonstrate the finiteness of the summation over n. Expression (|22p should be, by construction, finite. However, as 
one can immediately notice, the contribution of each term in Eq. (|22p is not separately finite. We need to combine 
them before taking the summation so as to cancel with each other. In the following we will demonstrate the finiteness 
of (|22[) explicitly. This, in principle, is not necessary, since the point-splitting regularization with subtraction of the 
appropriate counter-terms guarantees the finiteness of the renormalized Green function. However, aside from being 
a non-trivial check of the calculations, knowing explicitly how the various divergent pieces cancel with each other 
suggests a convenient strategy for the subsequent computation. 

First of all, let us rearrange the renormalized Green function by adding and subtracting its WKB counterpart: 



where 



\X))=G { r\x,X) 



k 2 a 



Stt 



2 {-Ai-Aa + To + Si + Sa} , 



(23) 



E x = £(i + l/2) 



1=0 



x 3 W (l) x 3 W (l) 



2EEC + 1/2) 



6 ) ^ l w„ 



n=l 1=0 
1 



X 3 W n {l) X 3 W„(l) 



(24) 
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and 



T 



E 

z=o 



Z + 1/2 1 



X 3 W„(/) Xy/J 



(25) 



The terms other than £2 are manifestly finite. The term £2 is cumbersome to evaluate, although straightforward. 
Once again, we can proceed in the standard way. Making use of the Abel-Plana summation formula, we rewrite T n 
as 



X 3 



1 1 



4W„(0) 



z + 1/2 
W n {z) 



-7j) dz - 



x 2 \ 



+ i 



dz 



+ 1/2 -iz + l/2 s 



1 V W n (iz) W n (-iz) 



(26) 



The contribution to £2 from the first term in the 
square brackets of Eq. l|2"6"|) is 



1 °° 1 

Pi = — Y — 



00 

Ot-3 E 



■2.r' ^ W n (0) 

jKg) \ 

— ^ |2Zi + (acoi - 2a 2 - 3a 3 ) Z 3 
[~^--2a 2 j (2a + xa') 



-03 3tT 



X4 " 2 1 >\ 2 -7 

■ [2a + xa ) Z-j 



(27) 



where we have substituted 1/W n (z), given by Eq. (|T4"|) . in 
the first equality and the explicit expression of 'J, given 
by Eq. (| 1 5[> . in the second equality. Here we have intro- 
duced generalized ^-functions of Epstein-Hurwitz type, 
defined by 



00 

E 

n=l 



(ujI+xMx))-^, 



(28) 



which occur in the computation of functional determi- 
nants of Laplacians on Riemann spheres (Some explicit 
applications can be found, for instance, in Refs. [40ll4ll].) 

We refer to the n-th term of the divergent piece as 
div[P] n . To be precise, here we define div[7- , ] n by the 
two leading terms in the large n expansion in the form 
of b\LJ n + 62 w" 1 , where 61 and b 2 are constants. In the 
above expression (|27[) only the term with Z\ is divergent. 
Hence, we can easily extract the divergent part of V\ as 



1 



2x 2 ui n 



(29) 



All the sub-leading terms are finite. 
The next term gives 



tE 



n=l 



1/2 



W n (z) 77 



dz 



2V7 



= P 2 ,l+P2,2, 

where we have defined 



(30) 



n=l 

00 -oo 



z + 1/2 x 2 



dz 



2V7 



Pi 



.7^3 j 



(z + 1/2)^ 

$3/2 



dz 



Here the integrations over z can be performed easily. Af- 
ter integration over z, P 2 ,i reduces to 



P2.1 



/ 



Z-i. 



(31) 



Hence, we can extract the divergent part of V2.1 as 



2 ( a 
div[P 2 ,x] n = - 7 ^ n - — 



(32) 



Also, we perform the integration over z in 7^2,2, substi- 
tuting the explicit expression of 4 1 , to obtain 



-x 3 (2cr + a;cr')Z 3 j 



02 f, / „„ 18a;/' x 2 f' 2 \ „ 



30/ r 



/ 

*/' 



,f 2 



Ax 2 [ 7 - -j- ) (2a + xa') Z 3 



+3a; 4 {2a + xa') Z 5 



"3 



(6x 2 <7 - zV) Z 3 



(33) 



Since the divergence in the above expression is contained 
only in the terms proportional to Z±, we can easily ex- 
tract the divergent part as 



1 f 13/ , 5f 

x 



(34) 
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where we have substituted explicit forms of a's given in 
Eqs. (HI and (S). 

The last term to evaluate is 



.r 



■ oo r .oo 

I 

3 



n=l 



dz 



iz +1/2 ~iz + 1/2 



' -'^/o e 2 ^-lV^„(iz) W n (-iz) 



Due to the exponential fall-off, the dominant contribution 
to the integral comes from the z ~ region. Hence, we 
can evaluate this term by expanding the part enclosed by 
the parentheses as 



iz + 1/2 -iz + 1/2 
W n (iz) ~ W n (-iz) 



oo 

= —i C n jZ 

i=i 



2i-i 



(36) 



The convergence of this series will be fast for large n, 
while it will be slow for small n. Therefore we divide Vs 
into two parts as 



^3 = ^3,1+^3,2, 



(37) 



where 



2 

n=l 



dz 



3 / „2ttz 1 

G i=l 



2i 



3,2 



E 



is + 1/2 

^ J e 2 ^ - 1 } V W n (iz) 



-iz + 1/2 
W„(-iz) 



E 



,(38) 



with an appropriately chosen value of j ma x- Then, on 
one hand, Vz,i can be evaluated numerically. As long as 
jmax is sufficiently large, the summation over n converges 
rapidly. On the other hand, we can perform the integra- 
tion over z in Vs i analytically term by term using the 
formula^ 00 z^^e 2 ™ - lj^tfe = T(2j)( R i2j)/{2n) 2 3 = 
(—iy~ 1 B2j/4:j, where £r is the Riemann ^-function and 
Sj = 1/2, B 2 = 1/6, B 3 = 0, B 4 = 1/30, B 5 = 0, B 6 = 
1/42, i?7 = 0, i?8 = 1/30, • • • are the Bernoulli numbers. 
Finally, we obtain 



•EE 

n=l j=l 



4j 



2,y 



(39) 



Since c„j is (^(l/w 2 - 7 " 1 ) for large n, only the part with 
j = 1 is divergent. Thus, using c n \ ~ —2/W n (0) and 
B2 = 1/6, we have 



div[7> 3 L = - 



6x 2 o 



(40) 



It is now a matter of trivial algebra to combine together 
(HHJ), d22), (O and ([30"]), to show that they cancel the 
contribution from the last line in l]22p. leaving a well- 
behaved expression. 

The same procedure applies in computing To, too. We 
just need to eliminate the summation over n setting n to 



zero in the expressions for Vs. For the evaluation of the 
counter part of V3 in computing To, we do not need to 
divide it into two pieces like V-s,\ and Vz,2 since there is 
no infinite summation over n. 

The remaining task is composed of two parts; numeri- 
cal calculation and analytic summation over n. The for- 
mer is necessary for A2, To, but infinite summation 
over n is unnecessary for these terms. Only the leading 
finite number of terms give sufficiently precise approxi- 
mation. The other terms, i.e. Ai and the other 7-"s, can 
be handled completely analytically. Moreover, from the 
above demonstration, we find that all the summations 
over n which appear in Vs are written in terms of Z n 
and its counter-terms in the form y\ b±u) n + 62^ 1 ■ <Z_i 
and Z\ contain divergences, and hence they require the 
counter-terms. The coefficients of the counter-terms, 61 
and 62, are appropriately chosen so as to cancel these 
divergences. Hence, we can simply define a regularized 
Z-i and Z\ by including the counter-terms proportional 
to X) n w " or En^' so tna * t ne y are nn ite. We denote 
such regularized Z_i and Z\ as 2+i and Z%, respectively. 
Then, summing up the regularized pieces gives the cor- 
rectly renormalizcd value. 



B. Summation over n 

In the previous subsections, we performed the sum- 
mation over I in the renormalized expression for ((f) 2 ), 
Eq. (f2"3"]) . The expression is composed of two parts; one 
is the part that requires numerical evaluation and the 
other is the part that requires infinite summation over 
n. In this subsection we will carry out the evaluation of 
these expressions. 

It is straightforward to perform the necessary numer- 
ical computations. The term Si requires evaluating W n 
and hence the exact mode functions, pf(x) and qf{x). 
For this term, summations over I and n are also neces- 
sary. In the practical numerical computation we need to 
truncate these summations at finite I and n. However, 
this term consists of the difference between W n and its 
WKB approximant W n , and hence it is the remainder of 
the approximation in this sense. Since the WKB approx- 
imation becomes better for large I and n, the convergence 
is basically fast. In the present approximation truncated 
at the next-to-leading order in the WKB expansion, the 
remaining terms cause the error in the estimation of (<f) 2 } 
of order 0(l~ 2 K , n m 2 x ). As we increase the order of WKB 
approximation, the error decreases more rapidly. How- 
ever, the price to pay for a faster numerical convergence is 
that the computation of the analytic part becomes more 
complicated. For small x, the above argument is com- 
pletely correct. For larger x, however, is dominated 
not by w but by a in wide range of I and n. When a 
is dominant, WKB expansion is not a good approxima- 
tion at all. (WKB series does not converge in this case.) 
Therefore even if we use higher order WKB approximant, 
we need to sum the difference up to very large I and n. 
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In this case, in order to reduce the computational cost, 
we can use the interpolation of the values of difference 
evaluated only at sparse grid points of I and n. Vz,i also 
has a summation over n but again it is composed of the 
difference between the expression (1361) and its truncation 
at j = jmax- As in the case of Ai, as we increase j max , 
the convergence of the summation over n becomes faster. 
Due to this property of fast convergence, no special trick 
is necessary in numerical evaluation of this term. 

The latter part should be treated fully analytically due 
to the presence of infinite summation. However, it is 
not so difficult because they arc all already written in 
terms of generalized zeta- functions, as we have seen in the 
previous subsection. These generalized zeta-functions, 
Z q , are functions of 



a" 

We need to evaluate them as a function of v 2 only once at 
the beginning of the whole calculation. Hence, the direct 
sum is one possible way to compute them. However, there 
are slightly cleverer ways to reduce the computational 
cost, which we will explain below. 

One way to truncate this summation is as fol- 
lows. Formally we write the summation defining Z q as 
S^Li F q {n). This summation can be directly performed 
numerically up to a certain value n*. For the remaining 
part of summation, we simply expand F q (n) in terms of 



-2s 



F,(n) = (ow)-«X;^fc© 



2fe+l 



(41) 



k=l 



where the coefficients d q k are independent of a, v and n. 
Then, we have 



oo n* 
n— 1 n— 1 

fcmax 

+ o-' £ d qk ((2k + l)v 2k -" +1 , (42) 



fc=i 



where 



E J j = cw-E^- 

z — ' n s z — ' n s 

n— n max + l n— 1 

In practice, we need to truncate the summation over k at 
a certain value fc max , which is not necessarily so large as 
long as n* is sufficiently large. In order that the expan- 
sion (|4"Tj) converges for v n > n„, n 2 must be larger than 



Hence, this method works efficiently unless 



is 



very large. 

On the other hand, for a very large positive value of v 2 , 
we can make use of a simplified version of the Chowla- 
Selberg formula (see for example [iol. |42|). 



S(s,p) 



(n 2 



2tt s 



Vi r( S -l/2) !_ 2s 

2 r( s ) 9 



(43) 



T(s) 



-s+l/2 



J2p S ' 1/2 K s - 1/2 (2npp) 



P =i 



The above identity basically rearranges the original sum 
in a convenient way, as the summation on the R.H.S. 
of the above equation converges very rapidly due to the 
exponential fall-off of the modified Bessel functions unless 
p is extremely small. 

Among Z q , Z^\ requires the most delicate treatment. 
By definition, 



Z_i = a f \J n 2 + v 2 — n — — J 



(44) 



The Chowla-Selberg formula does not directly apply to 
this expression since the R.H.S. of (|43|) diverges for 
s = 1/2. Of course, it is due to this divergence why 
we subtracted counter terms. The trick is to regularize 
the above expression as 

ZlJ^ajr f(n 2 + v 2 y-~n 1 - 2 ^ (~ - e\ v 2 ^ 1 ^ 

n=l \ ^ ' 

= S(-l/2 + e,v 2 )-( R (-l + 2e) 

-(l- e ) A«(l+2e). 

This expression reduces to Z_i if we set e to 0, and it 
is finite for any non-negative e. Hence, the limit e — > 
gives Z-\. We can use the Chowla-Selberg formula, take 
the limit, and obtain an explicitly regular expression, 



Z-i 



1 
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l-2(In(t;/2)+7 B ) 



P =i 



if_i(27rpw)\ 



P 



(45) 



Rearranging the generalized zeta functions by means of 
the Chowla-Selberg formula, i.e. as a series of modified 
Bessel functions, proves to be efficient when the argument 
of the modified Bessel functions is real and not small. 
For iJi, the regularized expression is simply given by 

Zl cs = S(l/2 + e,v)-( R (l + 2e). 

Taking the limit e — > 0, we obtain 

..-if 1 



Zi = —a 



— + \u{v/2) + lE -2^2 K (2irpv) )(.46) 



P =i 



The other generalized zeta functions Z q with q > 1, 
which do not require further regularization, are evalu- 
ated directly applying the Chowla-Selberg formula. We 
find 



z q = er« 



2 2T{q/2) 
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r w 2 > p 2 Kq-i (ZTrpv) 



r(«/2) 



(47) 



Finally, we mention how to evaluate Ai. It can be 
evaluated in an analogous manner. The first and the 
second terms in have exactly the same forms as Z_i 
and Z\, apart from the replacement of v 2 with rh 2 f /a 2 . 
We write it here more explicitly as 



Ai 



2a ~ 
T 



RZi 



(48) 



v 2 — ►m 2 / /a 2 



VI. RESULTS 

So far, we have not specified the form of the metric. 
For the numerical evaluation, however, we need to fix the 
function f(r). To illustrate the results we will choose the 
case of a Sch-AdS black hole, with metric given by ([T]) 
and / by ([2]). For notational convenience we define the 
quantity t^h = 2Mm~ 2 k 1 which characterizes the size of 
the black hole. 

The numerical analysis does not present particular dif- 
ficulties and can be performed in a straightforward man- 
ner. The WKB part is expressed in terms of analytic 
functions and can thus be evaluated very easily. Slightly 
more involved is the computation of the sums. The WKB 
approximation to next-to-leading order ensures that the 
convergence goes relatively fast, as 0(l~ 5 ,n~ 5 ). There- 
fore, after obtaining the solutions for the modes numeri- 
cally, the summations can be computed directly. 

After performing standard numerical checks, the vari- 
ous terms contributing to (4> 2 ) can be combined and some 
illustrative curves are reported for some representative 
values in the case of conformally coupled scalars (£ = 1/6 
and m = 0) in Fig. [T] (For n = I = and very small val- 
ues of M, the function $ may become negative. In this 
case, the n = I = mode has to be treated separately in 
the single sum To). 

Due to the complexity of the calculation, it is instruc- 
tive to look at the behavior of the vacuum polarization 
far away from the black hole, where we can expect to 
reproduce the leading AdS space result, which can be 
calculated independently. In fact, in pure AdS space the 
Green function is given by 



G(X,X') 



2 
3^ 



3 — cosh 



(ky/d(X,X'] 



where d(X, X') is the geodesic distance between X and 
X' . By subtracting, as in the black hole case, the counter- 
terms and then taking the coincidence limit, one gets: 

<^S--^- (49) 

Our results for Sch-AdS case can be fitted in the asymp- 
totic limit as 



2 ) ~ Ci + C 2 /f ■ 



(50) 



with C\ = —jg^?- The first term represents the lead- 
ing contribution to the vacuum polarization. We have 
checked the universality of C\ to high accuracy. 

The coefficient C-2 of the second term is hard to an- 
ticipate before calculation. If we use the result for the 
vacuum polarization at a finite temperature in Minkowski 
space, setting the temperature to the local Hawking tem- 
perature ka/2n^yj, we obtain the estimate C% = 
However, our numerical result for Ci is much smaller 
than this value in general. Since the temperature at large 
distance from the black hole is very low, a typical energy 
scale for the excitation is much below the inverse cur- 
vature length. Thus, it is not surprising that the rough 
estimate based on the result in Minkowski space does not 
hold in the present case. 



VII. DISCUSSIONS AND CONCLUSIONS 

In this paper we developed a method to compute the 
renormalized expectation value of a quantum scalar field, 
with mass m and coupling to the curvature £, in a ther- 
mal state on a spherically symmetric, asymptotically AdS 
black hole geometry. 

We followed the approach of Refs. 0,[H[j and employed 
the analytic WKB approximation and the point-splitting 
regularization to construct a regular expression for the 
coincidence limit of the Green function. We explicitly 
demonstrated the regularity of the Green function, and 
this allowed us to perform the renormalization term by 
term. 

Analogously to the asymptotically flat case, the WKB 
approximation arranges the vacuum polarization in a 
WKB-part' plus a remainder. One term depends on the 
WKB approximants and it can be evaluated analytically, 
although the computations become more cumbersome as 
the order of the approximation increases. 

It is very interesting to notice that the WKB approxi- 
mation organizes the analytical part of the vacuum polar- 
ization as a series of analytic functions. These functions 
take the form of generalized zeta functions, that occur in 
the computation of functional determinants of Laplacians 
on Riemann spheres. The coefficients of the expansions 
can be calculated order by order in the WKB expansion. 
We explicitly showed this to next-to-leading order in the 
WKB approximations. The other term is a remainder of 
the WKB approximation, in the sense that it depends on 
the difference between the exact solutions for the modes 
and their WKB counterpart. 

In the end, both terms have to be handled numeri- 
cally. The generalized zeta functions can be easily eval- 
uated using direct summation when their argument is 
small. When the argument is not small, it is convenient 
to rearrange these functions by means of a simplified ver- 
sion of the Chowla-Selberg formula [40| . The rearranged 
expression contains an infinite series of modified Bessel 
functions, but it converges exponentially fast when the 
argument is not small. The 'remainder' can be evalu- 
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FIG. 1: The figure illustrates the behavior of ((f> 2 ) for conformally coupled (£ = 1/6 and m = 0) fields. The curves refer to: 
M = 100, 50, 5, 0.5. 



ated by directly solving the mode equation numerically 
and taking the difference with the WKB counterparts. 
These terms converge very rapidly, as we expected from 
the fact that the WKB approximation works very well: 
at next-to-leading order the remainder is 0(l~ 5 ,n~ 5 ). 

In the last part, we explicitly illustrated the method by 
computing the vacuum polarization for conformal fields 
on a Sch-AdS black hole background. We finally dis- 
cussed the asymptotic behavior of the vacuum polariza- 
tion and compared it with the pure AdS result which 
can be calculated independently. We find that our result 
reproduces the leading universal (constant) behavior to 
high accuracy. 

We are now using the methodology presented in this 
paper to compute the energy momentum tensor for a 
quantum field on an asymptotically AdS black hole ge- 



ometry and we hope to report on this soon. 
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